Electron localization in one dimension obtained from combined exact 

diagonalization - ab initio approach 



O 
O 

(N 

G 

in 



■4— » 

o3 



C 

o 
o 



> 

On 
OO 
O 

\o 

O 



a 

I 

C 

O 

o 



X 



Jozef Spalek and Adam Rycerz 
Marian Smoluchowski Institute of Physics, Jagiellonian University, ulica Reymonta 4, 

30-059 Krakow, Poland 
(June 4, 2001) 

Exact ground-state properties are presented by combining the diagonalization in the Fock space 
(and taking all hopping integrals and all two-site interactions) with the ab initio optimization of the 
Wannier functions. Electrons are essentially localized for the interatomic distance R ~ 2.0 A for 
s-like states, when the quasiparticle mass is divergent. The momentum distribution dispersion is 
proposed to define the localization order parameter. Dimerization and zero-point energies are also 
discussed. The method provides convergent results for N > 8 atoms. 



PACS Nos. 71.10.Fd, 71.15.Fv, 31.25.Nj 

One dimensional systems range from organic metals |l]] 
to quantum rings and wires [|| , and to nanotubes . In 
their description the role of the long-range Coulomb in- 
teraction is crucial because of reduced dimensionality, for 
which the charge screening becomes less effective Q . The 
existing exact solutions of the parametrized models with 
inclusion of intersite interactions prove the existence 
of the metal-insulator transition for the half-fillcd-band 
case, in contradistinction to the corresponding Hubbard- 
model solution Q , for which the system is insulating even 
for an arbitrarily small Coulomb repulsion. The exis- 
tence of such metal-insulator the transition has been also 
discussed within the density-matrix renormalization 
group (DMRG) method when the second-neighbor hop- 
ping is included. A separate question concerns the ap- 
pearance of the Tomonaga-Luttinger behavior || in the 
metallic state, for which some evidence has been gathered 
H . In brief, the search for a proper description of those 
systems as quantum liquids (and their instabilities) is one 
of the basic problems in the physics of low-dimensional 
systems. 

In the above theoretical analysis [1-8] the solutions 
have been disscussed as a function of the microscopic pa- 
rameters, which are not easy to measure. Therefore, one 
assumes that they should be determined first from a sep- 
arate single-particle approach. In following this route 
one must avoid counting twice the interaction, as dis- 
cussed carefully in the papers implementing the LDA+U 
and SIC Q methods. We have proposed fjj] a 
new method in which the single-particle (Wannier) wave 
functions are allowed to relax in the correlated state and 
thus are determined by optimizing the exact ground state 
energy obtained from the diagonalization in the Fock 
space. In other words, we treat properly the interac- 
tions first due to their strongly nonperturbative nature 
and then readjust the single-particle functions {wi(r)} 
by setting Euler equation for them, which plays the role 
of the renormalized wave equation for a particle in the 
correlated ground state. In such approach the problem 



of double counting the interaction does not arise at all. 
Additionally, we include all the hopping integrals and all 
two-site interactions to make the solution more complete. 
The method is executable on a desktop server for the 
number of atoms N < 12. What is probably the most 
remarkable formal feature of these calculations is that 
the electronic correlations make the wave function more 
tightly bound to the atoms and thus limiting the interac- 
tion range. In effect, the results are converging very fast 
for N > 8 atoms meaning that the optimized Wannier 
functions and the interaction parameters are relevant at 
most up to the third neighbors. The results are detailed 
below and to the best of our knowledge they represent 
the first analysis of the strong correlation effects as a 
function of the lattice parameter within the exact acount 
of both the interaction and the wave function without 
limiting the range of either the hopping processes or the 
two-site interactions. The obtained results show that the 
method is particularly useful for accurate studies of elec- 
tronic correlations in quantum dots and rings. 
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FIG. 1. Momentum distribution n^a for electrons in linear 
ring of N — 10 atoms; the interatomic distance R is specified 
in units of Bohr radius ao . The continuous line represents the 
parabolic interpolation, which is of the same type for both 
k > Icf and k < ftp. 
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The nature of the electron momentum distribution |13|] 
rifc ff = < a\ a a,ka > is shown in Fig. 1 for N = 10 atoms 
(we use the periodic boundary conditions). For the in- 
teratomic distance R — 3ao, where ao is the Is Bohr 
radius (setting the length scale), this is essentially the 
Fermi-Dirac function with a tail extending to the Bril- 
louin zone boundary. The smearing out of the distribu- 
tion with increasing R suggests that the electronic states 
transform from itinerant to localized states, as exempli- 
fied for R — R c 3.929. The continuous line represents 
the parabolic fit of the same type for both k < kp and 
k > kp. A natural question to be dealt with is whether 
the Fermi ridge disappearance at R = R c is reflected in 
any singularity in other properties. 
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FIG. 2. Correlations function versus distance R, depicting 
the crossover from itinerant to localized state (see main text), 
for N = 6 to 10 atoms. The shaded areas are drawn to illus- 
trate the convergence of the results in the large R limit. 

To distinguish quantitatively between those states we 
plot in Fig. 2, the following basic quantities (as a func- 
tion of R) defined in both itinerant and localized states: 
(i) the site spin magnitude 0m = (4/3) (S?\, where 

Si = (Sf,Sr,Sf) = (o] T o»;, a^ajf, (n it -n a )/2) 
is the electron spin on site i, (ii) the spin correlation 
function Oaf = - (S, • S i+ i), and (iii) 9 M i = ^o 2 {n k <r}, 
where a 2 {n^a} is the dispersion of the statistical distri- 
bution defined as: 




The averages are for the ground state, which is deter- 
mined via an exact diagonalization in the Fock space. 
The quantity (9m = 1 — 2 < >) takes the value 

(1/2) in the ideal gas limit and approaches unity in the 
atomic limit, where we have a Pauli spin on each atom. 
Oaf approaches the value (3/4) for the singlet config- 
uration of atomic spins, whereas a 2 {nfc CT } acquires the 
value 1 in the gas limit (rifc CT = O (fj, — e&)) and van- 
ishes for an even momentum distribution (nfc CT = 1/2), 
when the particle position is sharply defined on atom, i.e. 
for the localized electron states characterized by atomic 
states. Thus, the quantity Omj plays the role of the order 
parameter for this crossover behavior, since it clearly dis- 
tinguishes between the complementary momentum and 
position representations of the system quantum states. 
From Fig. 2 it follows that for R/ao = 5 all three pa- 
rameters acquire (within 5%) their asymptotic values for 
purely atomic states. 
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FIG. 3. The bandwith-to-interaction ratio versus R; the 
Hubbard point for localization [15] is marked. The inset pro- 
vides the optimal size a -1 (in units of ao for s atomic orbitals 
composing the optimized Wannier function. 

The analysis of metal-insulator transition induced by 
electron-electron interaction started with the works of 
Mott JIJ] and Hubbard |J|. The first of them intro- 
duced the critical carrier concentration for the transition 
to take place in a discontinuous manner, whereas the 
second introduced the critical bandwidth to interaction 
ratio (of the order of unity). Our method of approach 
can be used to relate the above criteria, as we discuss 
next. Namely, to relate our method to the original ideas 
of Mott and Hubbard, we have plotted in Fig. 3 the ratio 
of bandwidth W = 2| Y^jU) *y l> wnere Uj is the hopping 
integral between the neighbors (ij) (calculated through 
the relation tij = (wi\H\ \wj) for the basis {wi} opti- 
mized in the correlated state and for full single-particle 
potential in Hi), to the effective short range Coulomb 
interaction U — K, where U = (wiWi\ V12 \wiWi) and 
K = (wiW i+ i\V l2 \wiW i+1 ). The value W/ {U - K) = 1 



2 



marked in this figure represents roughly the dividing line 
between metallic and Mott insulating states for three- 
dimensional systems [ p"5| . This point, achieved for R « 
2.7ao, does not reflect any characteristic point for our 
system. Instead, the localization is practically achieved 
for the distance about twice as large, as is the size of the 
atomic states composing the Wannier function [Q (char- 
acterized by the quantity a , see the inset), which nears 
its atomic value for the s-like state. Again, the results 
for N — 8 and N = 10 are very close to each other; this 
is the reason why we have shadowed the areas between 
the corresponding curves in both Figs. 2 and 3. 

On the basis of Figs. 2 and 3 we can estimate the lo- 
calization threshold for our system. The corresponding 
Mott criterion Jl4j , generalized to d dimensions takes the 
form n~J d aH ~ 0.2, where n c is the carrier concentra- 
tion, and an is the size of the states at the localization 
threshold. In our situation of neutral chain: n c = 1/R, 
and an = oT 1 so that for R w 5ao this criterion takes 
the form (ao/R) (a _1 /a ) « 0.95/5 « 0.2, a suprisingly 
close value to that of Mott (which reflects the long-range 
nature of the Coulomb interaction). Thus, our results 
provide a support for the Mott criterion rather than for 
the Hubbard one [fl5|| . In other words, the metallicity 
extends well beyond the W = U limit and this must 
be caused by the inclusion of more distant hopping pro- 
cesses. One may also say that the Mott-criterion univer- 
sality originates from the long-range nature of the inter- 
action, which imitates the higher-lattice dimensionality. 

Hence, the conclusion Q about the universality of the 
insulating state for the Hubbard chain does not extend 
to the Id models with a realistic account of the electronic 
structure. This conclusion is very important also because 
it removes one of the main objections against using the 
itinerant (or even effective mass) states for quantum dots 




LU 

a: 



100-I , , , , , , , , , , 1 

3,0 3,2 3,4 3,6 3,8 4,0 

R/a 

o 

FIG. 4. The critical behavior of the quasiparticle mass at 
the Fermi level: The computed points are fitted with the curve 
discussed in the text. 

We can determine directly the effective mass at the 



Fermi level. Namely, the calculated band mass uif is 
about 40% enhanced near the localization threshold. The 
quasiparticle mass m* F is found from the relation m F = 
(Anfc F ) _1 mi?, where the first factor is the usual Fermi- 
liquid enhancement (the discontinuity at the Fermi 
level marked in Fig.l). The results are shown in Fig. 4. 
A clear critical behavior is detected: m* F = A\R — i? c |~ 7 , 
with A = 10.2, R c = 3.92, and 7 = 4/3. This quan- 
tum critical behavior is obtained, since we emulate the 
discrete distribution rik with a continuous parabolic in- 
terpolation when determining An& F . The localization 
threshold R c is about 10% higher for N = 8. It is tempt- 
ing to speculate that with the increasing N — > 00 the 
Mott and the Hubbard criteria for R c may coalesce. The 
masses near R = R c are huge and reach those known 
only for heavy- fermion compounds. 

The one-dimensional systems are unstable with respect 
to the dimerization jjj. We have determined the ground 
state energy of such a state on the same level of precision 
as above and have additionally minimized it with respect 
to the amplitude of the lattice distortion. In Fig. 5 (bot- 
tom) we have compared the energy contribution AE due 
to the dimerization with that due to the zero-point mo- 
tion (top) . Both the dimerization amplitude and AE are 
strongly reduced for R/ao > 5, a feature inevitably con- 
nected with the long-range nature of Coulomb interac- 
tions, which drives the system towards spatially periodic 
state (the Mott insulating state on the lattice replaces the 
Wigner-crystal state S for the electron gas). Therefore, 
our results and conclusions above should remain intact 
for longer chains. Additionally, the zero-point energy 
overcomes the dimerization energy for a light (hydrogen) 
chain and has been estimated in the following manner. 

In the harmonic approximation, the phonons have en- 
1/2 

ergy u>k — 2 (C/M) 1 sin (irk/N), where M is ionic mass, 
and the elastic constant is calculated by a numerical dif- 
ferentiation C — N~ 1 d 2 E G /dR 2 . Formally, these for- 
mulae are valid only if the energy has an absolute mini- 
mum. Here we assume that the system is closed in a box 
of length NR and thus the global repulsive interaction 
between the atoms (the situation with dEc/dR < for- 
given R) is balanced out by the environment. Also, as 
the k=0 mode is a Goldstone mode, we include only the 
modes with fc 7^ 0. The phononic contribution to the 
ground state energy (in the atomic units) is then 

ae ° ={m) bwj £ sn w (2) 

where m is the bare electron mass. On the basis of the 
relation for the fc-th mode (1/2) Muk (ARk) 2 = 
where [ARk) is the mode contribution to the zero-point 
vibrations, we can easily estimate the total amplitude 
(AR) 2 , which in the atomic units has the form 
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In the A 7 — > oo this result give (Ai?) ~ log A 7 ", providing 
the dynamical lattice instability in one dimension. At the 
localization threshold (i.e. for R/clq pa 5) and for N = 8, 
we have that Ai? « 0.12, a substantial. By comparison, 
the shift due to the dimerization is « 0.06 . Also, the 
Peierls distorted state extends to the localized state (up 
to R pa 6.5 do). 
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FIG. 5. The zero-point (top curves) and the dimerization 
(bottom curves) energies for N — 4-^10 atoms versus R (a~ in 
is the optimal size of the atomic orbitals) . 

In summary, we have determined the microscopic cri- 
teria for the transition from the itinerant to the localized 
states in a one-dimensional system of a finite size and 
have illustrated those findings on the example of a corre- 
lated quantum ring. The new method of optimizing the 
single-particle wave functions in the correlated state fl^ l 
proves thus valuable in the exact treatment of nanoscopic 
systems. We have calculated all the properties as a func- 
tion of the lattice parameter. The work was supported 
by KBN Grant No. 2P03B 092 18. 
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